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In the pre-reheating era, following cosmic inflation and preceding radiation domination, the energy 
density may be dominated by an oscillating massive scalar condensate, such as is the case for 
V = m 2 (j> 2 /2 chaotic inflation. We have found in a previous paper that during this period, a 
wide range of sub-Hubble scale perturbations are subject to a preheating instability, leading to the 
growth of density perturbations ultimately collapsing to form non-linear structures. We compute 
here the gravitational wave signal due to these structures in the linear limit and present estimates for 
emission in the non-linear limit due to various effects: the collapse of halos, the tidal interactions, the 
evaporation during the conversion of the inflaton condensate into radiation and finally the ensuing 
turbulent cascades. The gravitational wave signal could be rather large and potentially testable by 
future detectors. 



I. INTRODUCTION 

One of the most salient features of inflationary cos- 
mology is undoubtedly the generation of both scalar and 
tensor modes of metric perturbations on super-horizon 
scales through the amplification of vacuum fluctuations. 
Scalar modes correspond to genuine density fluctuations 
which are probed with ever increasing accuracy by the 
measurements of the temperature fluctuations of the cos- 
mic microwave background and large scale surveys. Ten- 
sor modes, however remain to be detected. Current (e.g. 
Planck [H) or next generation experiments (e.g. BBO (2|, 
DECIGO 3) should in principle detect the relic gravita- 
tional wave background if the energy scale of inflation is 
close to the GUT scale. As is well known, gravitational 
waves are much harder to detect because of the intrin- 
sic weakness of gravity, but at the same time and for 
the same reason they offer an invaluable probe of early 
Universe physics. 

In the past two decades, it has been realized that grav- 
itational waves could also be produced in the primordial 
Universe and on sub-horizon scales due to a rapid distur- 
bance of an otherwise equilibrium homogeneous state. 
This may occur for instance in first order primordial 
phase transitions 043 1 through bubble collisions and the 
dissipation of energy in turbulent cascades or 
via the interactions of waves generated by parametric 
amplification in preheating scenarios involving multiple 
scalar fields or tachyonic couplings of the inflaton [13l - [l8j . 



'Electronic address: jedamzik@lpta.univ-montp2.fr 
t Electronic address: lemoineC 
t Electronic address: jmartinC 



These various scenarios lead to different predictions at 
different frequencies depending on the detailed processes 
at work, with interesting prospects for detection [l9j |. 

In a previous paper (hereafter Paper I) [20j, we have 
shown that single field inflation generically leads to the 
growth of density fluctuations in the pre-reheating era for 
a certain range of wave numbers on sub-Hubble scales, 
and we have argued that one of the possible signatures of 
this effect would be the generation of gravitational waves 
(GW). The objective of the present paper is to calculate 
explicitly the spectrum of gravitational waves emitted by 
the interactions of these density perturbations. In order 
to avoid unnecessary model dependencies, we focus on 
the simplest model of inflation, namely chaotic inflation 
with potential V(tf>) — m 2 (f> 2 /2. Then, the normalization 
of the density perturbation power spectrum to the value 
measured by cosmic microwave probes sets the scale for 
the mass m and the Hubble parameter at the end of infla- 
tion. The only parameter left in the model is the reheat- 
ing temperature of the Universe, T r h, which marks the 
beginning of the radiation dominated era and which is 
directly related to the decay constant of the inflaton, 
T r h ~ 0.2-y/r^A/pi (with Mp\ the reduced Planck mass). 
In the following, we calculate the present day gravita- 
tional wave spectrum in the linear regime and provide 
detailed estimates for gravitational wave emission due to 
non-linear effects, which are associated with the growth 
of the density contrast on different scales beyond unity. 
We find potentially observable signals, all the more so 
when the reheating temperature is lowered. A lower re- 
heating temperature is associated with a longer epoch 
during which density fluctuations may grow as well as a 
frequency closer to the range of maximum sensitivity of 
next generation gravitational wave detectors. 

Recently, the problem of gravitational wave emission 
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by the growth of density fluctuations in an early matter 
era has been discussed by Assadullahi and Wands plj ]. 
The problem at hand is very similar but there are sev- 
eral differences between this study and the present one. 
Firstly, Assadullahi and Wands do not specify the mech- 
anism through which fluctuations can grow but assume 
that all sub-horizon modes are unstable if the equation of 
state of the background field is on average pressure-less, 
as happens for m 2 <f) 2 during the inflaton oscillatory phase. 
As demonstrated in Paper I however, growth of fluctu- 
ations is a preheating effect, which is limited to a band 
of wave numbers, whose range depends on the precise 
shape of the potential. Borrowing from this calculation, 
we are thus able in the present paper to relate the growth 
of fluctuations with the initial post-inflationary power 
spectrum of metric fluctuations. More importantly, we 
provide in the present paper detailed estimates for the 
non-linear regime of matter fluctuations, while the au- 
thors of Ref. [2l| have focused their study only on the 
linear regime. 

The outline of the paper is as follows: In Section |TT] 
we recall the main findings of Paper I and set the stage 
for the calculation of the gravitational wave signal done 
in Section HTT1 in the linear limit. Section llVl discusses in 
four subsections, gravitational wave emission in the non- 
linear regime, due to the first collapse phase of density 
fluctuations, due to the non-linear phases after first col- 
lapse, due to the evaporation of the collapsed halos at 
reheating and due to the dissipation of energy in turbu- 
lent cascades after reheating, respectively. Conclusions 
are drawn in Section W\ 



II. GRAVITATIONAL INSTABILITY IN THE 
PRE-REHEATING ERA 

In this section, we briefly recall the main findings of Pa- 
per I. In that article, we have discussed the evolution of 
the Mukhanov variable t>k and of the curvature perturba- 
tion £k during the pre-heating epoch. These two variables 
control the dynamics of the scalar metric perturbations. 
During pre-heating, the inflaton field oscillates according 
to 



/ ^end \ "V^ 

<f>(t) ~ 0cnd ("^"J sin (mi) , 



(1) 



where cn( j — mpi/(2v3) is the inflaton vacuum expec- 
tation value at the end of inflation, and where mpi = 
s/SttMpi. The mass m is fixed by the COBE/WMAP 
normalization, m ~ 1.7 x 10 13 GeV. The Mukhanov vari- 
able Vk is related to the curvature perturbation by 



(2) 



where t\ = —H/H 2 is the first slow-roll parameter ex- 
pressed in terms of the Hubble parameter H — a/ a. The 



quantity a denotes the Friedman-Lemaitre-Robertson- 
Walker scale factor, k = Mp 2 and a dot means a deriva- 
tive with respect to cosmic time. 

In Paper I, we have shown that there exists an instabil- 
ity regime of narrow resonance provided the wave number 
of the Fourier mode is such that 



< - < y/SHm . 
a 



(3) 



Once a mode has entered this resonance band, it remains 
inside it during the whole oscillatory phase. Moreover, 
we have also demonstrated that, in the resonance band, 
Vk oc a and, therefore, given the definition of this quan- 
tity given above, that remains constant. Note that 
usually the curvature perturbation remains constant only 
on super Hubble scales while decaying on sub Hubble 
scales, unless there is gravitational instability. By using 
the perturbed Einstein equation 

k 2 



Sk = —- 
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with the shorthand notation ^(77) = Spk/p introduc- 
ing the fractional mass-energy density perturbation, and 
given that p oc I /a 3 , we may see that Sk grows as a(t) 
on sub-Hubble scales and remains constant on super- 
Hubble scales. The preheati ng ins tability for all sub- 
Hubblc modes with k/a < v / 3i?m may be understood 
as the gravitational instability of the oscillating scalar 
condensate. 

As already stated, the amplification described above 
occurs only for a limited range of scales, k £ [k m - ln , /c ma x] 
in terms of comoving wave number. We now discuss this 
point in more detail and specify fc m i n and fc max . A priori, 
the smallest unstable scale (in wavelength) is that which 
enters the resonance band at the end of reheating. How- 
ever, it has never been outside the Hubble scale during 
inflation and, therefore, the power stored on that scale 
at the end of inflation is suppressed. As a consequence, 
for practical purposes, one can consider that the small- 
est scale for which the amplification occurs is such that 
k ~ a cnd H oad , i.e. k = k/(a cnd H cud ) ~ 1, describing 
modes that leave and re-enter immediately the Hubble 
radius at the end of inflation, hence fc max = o, end H elild 
(fcmax = !)• More detail on this issue can also be found 
in paper I. 

On the other hand, the largest scale is the one for which 
modes re-enter the Hubble radius at the end of the mat- 
ter dominated pre- reheating epoch, i.e. k mm ~ a r h-ff r h, 
where H r h is the Hubble parameter at the end of the re- 
heating epoch or at the beginning of the radiation dom- 
inated phase. One can show that the number of e-folds 
between the Hubble radius exit of fc m in during inflation 
and the end of inflation is formally given by: 



In 



"•end 



= -*->-i 
2 2 



-p- Y k 2 



(•5) 



where W_i denotes the Lambert function. If one neglects 
the evolution of the Hubble scale over the interval be- 
tween Hubble radius exit of fc m ; n and the end of inflation, 
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the above can be approximated with afc min /a cn d — k min . 
In terms of the reheating temperature and inflationary 
scale given by 



rh 



end 



5.9 x 10 



5 ( T rh \ 

V 10 9 GeV J 



2/3 



(0) 



The last equation assumes g* = 230 degrees of free- 
dom at reheating. For a reheating temperature of T rn ~ 
10 9 GeV, one finds ln(a cn d/afc min ) ~ 11. Therefore, for 
T rn ~ 10 9 GeV, the scales we would be interested in in 
this article are those which left the Hubble radius during 
inflation between eleven and zero e-folds before the end 
of inflation. 

Finally, we also need the initial power stored on each 
scale k € [k min , /c max ] at the end of inflation. This is 
described by the power spectrum of the curvature fluc- 
tuation defined by 



V c (k) 
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2~n 



2 ICfcl 
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'pi 



(7) 



The corresponding primordial power spectrum for 
m 2 (/) 2 /2 inflation correctly normalized to COBE/WMAP 
observations is shown in Fig. [T] This figure has been ob- 
tained through a numerical integration of the evolution of 
Cfc (through Vk)- as discussed in Paper I, the slow-roll ap- 
proximation at first or even second order cannot reliably 
predict the shape of the power spectrum in the range of 
comoving wave numbers of interest to us, k < a C nd-ffcnd, 
which exit and re-enter the Hubble radius not too long 
afterward (in particular, before reheating is completed). 
In the notation of Fig. [U AT ex it — iV ond ~ — 10 — > 0, the 
lower bound depending on the exact value of the reheat- 
ing temperature. 



III. GENERATION OF GRAVITATIONAL 
WAVES AFTER INFLATION: THE LINEAR 
REGIME 

Let us now consider the production of GWs induced by 
the phenomenon described in the previous section. Our 
analysis follows the main thread of the analysis presented 
in Refs. [TJ1 [TJ| . We work with the metric 

ds 2 = a 2 {(l + 2$)dry 2 - [(l + 2$)5 ij +h ij ]dx i dx j } , (8) 

where <J> is the gauge-invariant potential (i.e. the poten- 
tial in longitudinal gauge) which describes the presence of 
inhomogeneities in the inflaton field and hij is the trans- 
verse traceless spatial tensor accounting for the gener- 
ation of GWs. We are interested in the amplitude hij 
of the transverse and traceless which obey the following 
equation (in Fourier space) 



k 2 - — 
a 



= 2naT^ 



(>7>*0 



(9) 



10" 



10" 



.* 10" 

Kr 

10" 



10" 



-60 



-50 



-40 



-30 
«(*)-N. 



-20 



-10 



FIG. 1: Primordial power spectrum of curvature perturba- 
tions for the large field inflation model with V(</>) = m 2 (j} 2 /2. 
Here N cxit (k) - iV cn d is In [a cross (k)/a cnd ], with a cross (fe) the 
scale factor at Hubble radius crossing of mode k during infla- 
tion and a on[ j the scale factor at the end of inflation. 



where we used the rescaled amplitude — ah^j and 
with Tij(r),k) = J dxc lk - x T lj {-q,x)/{2i:f/ 2 . The trans- 
verse traceless part of the spatial energy-momentum ten- 
sor T^ T is obtained by 



T'TT — I 
ij *-ij£m -L £m 5 



(10) 



with 



-Ly£m= Pii{ek)P 3 m{e-k) - 7;Pij(ek)Pem(e k ) , (11) 



and 



Pij (e fc ) 



(12) 



Here the quantity e k = k/k is the unit vector in the k 
direction. We are interested in the generation of GWs 
by the parametrically growing inflaton density pertur- 
bations between the end of inflation and the epoch of 
reheating, i.e. the epoch of inflaton decay. We will see 
below that the lowest non- vanishing contribution of these 
perturbations to T^ T appears only at second order in the 
perturbative quantity S k ■ We will thus work to this order 
in T^ T but only to first order in all other quantities. If 
we neglect the term a" /a, as appropriate for sub Hubble 
scales, then the solution of the above equation can be 
written as 



= Aij (k) sin [k (17 - r){)]+Bij (k) cos [k (77 - r?f)] 



hij(r), k) = 


Aij 


where 




Mj = 


2 k 


~~k 




2k 


Bij = 


T 



(13) 



drcos [k (m ~ t)] cl{t)T?? (t, k) ,(14) 



2k f m 

T JdT sin [k fa - t)] a(r)T^ (r, fc) . (15) 
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In these expression, we have assumed that the source is 
non negligible in the time interval r) ln i < r\ < r]{. In our 
calculation, the initial time is set by the time at which 
a given mode enters the Hubble radius, while r]f is given 



by the time at which a mode reaches non-linearity. A 
discussion of generation of GWs in the non-linear regime 
is deferred to the next section. The energy density of 
GWs is given by 



^£</^fo, *)%•(»>,*)>. (16) 



4«;a 4 (7y) 



where terms proportional to the Hubble constant have been neglected, appropriate for k/a 3> H. Using the Fourier 
expansion, one obtains 

*wfo) = 4^)£ / (^b / (^ e ^ (fe " fe ' ) ' a,fcfc '{(^( fe )^( fc ')> cos Mr,-*)] cos [k' (r,- m )] 
- <A y -(fc)B* (*')) cos [k fa - m )} sin [k 1 (ry - m )} - (B l3 (k)A* (k')) sin [k ( V - m)} cos [k' ( V - m )} 
+ (BaWBfjW)) sin [k (r, - nt)] sin [k 1 ( V - rj f )] j . (17) 

where we have also used Ajj(fc) = A|-(— fe) and Bij(k) = B*j(—k) imposed by hij(r],x) being real. Expression 
Eq. (TTT)) contains four terms. Let us first discuss the first one. One needs to express the correlator of two coefficients 
Aij. Using the expressions established before, see Eq. (Ti"4"|) . one is led to 

(A iJ (fc)^.(fc')) = — / dr dT'coB[k( m -r)}co S [k'( m -r')}a(T)a(T')(T^(r,k)T^(r',k')), (18) 

where one has 

(13* (r, k)T^*(r',k')) =Uj tm Uir. (T tm (T, fe)T; s (r', fe')) , (19) 

with the projectors defined in Eqs. (|11|) and (fT2|) . 

To go further, one needs to specify what the stress-energy tensor is. We take the stress-energy tensor of a pressure- 
less fluid since, when the field rapidly oscillates at the bottom of its potential, the pressure vanishes on average. 
Therefore, we have 

Tij = pa 2 ViVj , (20) 

where v represents the velocity of the particles in the pressure-less fluid. Note that Eq. (I2TJ1) is already second order 
in the perturbation since Vi is first order. The velocity can be calculated from the continuity equation which reads 

9 f g P\ ,iy« = 0. (21) 



9t V P 

Fourier transforming Eq. (1211) . and noting that Sp^/p oc a(t) for growing density fluctuations in which we are interested 
(as was discussed in the previous section), we may relate the Fourier amplitudes of density perturbation and velocity 
via 

v{ V ,k) = - l -^5 k {T 1 )e k . (22) 

Then, inserting the above expression for the velocity into the definition of the stress-energy tensor Fourier component, 
one obtains that 

(TUT,k)T r * s (r',k')) = ^p{r)a\r)H\r) P {r')a\r')H\ T ') J d 9 J ^f^Z^^^ 

x (6 g (T)6 k - g (T)6* p (T')6* k ,_ p (T')) . (23) 
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Then we use the fact that the density contrast grows proportional to a when the mode is inside the Hubble radius, while 
it remains constant from the end of inflation until Hubble radius crossing. Therefore, in the sub-horizon regime, one 
can write S q (r) = S q (r cnd )a(T) /a-t c , with ah c the scale factor at horizon crossing and 5 q (T end ) the value of S q at the end 
of inflation. In order to match both asymptotic behaviors in the super-horizon, namely 5 q (r) -> S q (T end ) = (6/5)Cq and 
in the sub-horizon regimes , 5 q {r) — > (2/5) [q/(aH)] 2 ( q , see Eq. ((4]), one defines ahc as that at which q = \/3ahc-ffhc- 
Given that aH oc a -1 / 2 in this preheating stage, this indeed implies that [q/(aH)} — 3(a/ah c ), hence the two 
asymptotes are matched at a = ai lc . Hence, one can write 

(Tt m (T,k)T rs (T',k')) = 7^X3/^32 J dq J d P1e(k m - q m )pr(K - p s ) (S q (T end )S k _ q (T end )5*(T end )Sl,_ p (T end )) , 

(24) 

The next step consists in evaluating the correlator. Using the Wick theorem in order to express the four-point 
correlation function in terms of three two-point correlations functions and discarding a homogeneous piece we are not 
interested in, one arrives at 

1 1 r 

(25) 



2 1 f 

(T em (T, k)T rs (T',k')) = J^j3j^j2 J dqq e {k m - q m )q r (k s - q s )<jq<j\ k _ ql S{k - k') 



where we have defined (5k(T cn d)5 k , (r cn d)) = 0fc<5(fc — fc'). Finally, we act with the projectors on the last expression 
and one obtains 

J2(W(T,k)T^*(r>,k')) = (2^3(3^2 / dgq 4 sin 4 aa 2 af fc _ q|( 5(fc-fc'), (26) 

where a is the angle between the vectors fc and q. Putting everything together, one can therefore evaluate the 
correlator between the coefficients Ay. We find 

a 2 rm rm 1 1 

^(^■(fc)^-(fc')) = drj^ dr'a(r) cos [fc ( m - r)] a(r') cos [fc' ( m - r')] ( ^ )3 ^ 

x Jdqq 4 S m A aa 2 q tf k _ ql 6(k-k'). (27) 
Therefore, the first term in Eq. (117[) is given by 

dfc »r dr r 

.2 

'q u \k-q\ 



K f dfc f f 

pii A) (V) = J 7^3 cos 2 [fc (77 - r)t)} I dr / dr'a(r) cos [fc ( m - r)] a(r') cos [fc (r) { - t')] 



1 1 

(2tt) 3 (3k) 5 



dq q 4 sin 4 a a 2 af k ,. (28) 



Taking the time average of the previous expression over a period replaces cos 2 [fc (rj — rjf)] with 1/2, and cancels the 

cross-terms (AB*), (BA*). The term p g w^ is similar to Pgw^ except that the term cos [fc (rji — r)] cos [fc' — r')] is 
replaced by sin [fc (rji — t)] sin [fc' (rjf — r')]. When combined, these two terms lead to a present day energy density 

k f f m f m 1 

^ (ro) = W^J dk lJ T J v J T ' cos[HT ~ T ' )]a ^ 

X J dq q 4 sin 4 a a 2 cx 2 fc _ q| = ^7^4 / dfc k 2 J(k)X(k) , (29) 

where the functions J and I are defined as follows 

/•Vi rVi r 
J(k) = / dr / dr' cos [k (r - r')] a(r)a(r') , I(fc) = / dqg 4 sin 4 a <^cf fe _ g| ■ (30) 

J 77i n i */ 77i n i «/ 



The first integral can be done exactly. The result reads 

, v 6 



= 77777 ( 77 j {8 + x 4 + , 



1111 



2 [4 + 4a; f x ini - 2x 2 ni + x 2 (x 2 ni - 2)] cos (x f - x ini ) + 4 (x ini - x f ) (2 + x f x ini ) sin (x { - a; ini )| , (31) 
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with x = 3t^jt 1/ ' 3 A;/aeiid = 2(a/a cn d) 1 / 2 fc/(a e nd-ffend)- Note that x > 1 for sub-Hubble modes. The second integral 
reads 

X(k) = 2tt r^dq [ + dA*a^af fc _ a | (l-^ 2 )V. (32) 



Using <5q(r en d) = (6/5)Cg and Eq. ([7]), we have 

36 2tt 2 



cr 



2 



9 25 g 



Pf(g), (33) 



where the power spectrum of the curvature perturbation is calculated at the end of inflation, see Fig. U As a 
consequence, one obtains 



= 8 * K*) JL * /- d " p<(,,p< ^* + 42 - ^^-^ ■ (34) 

Using the above expressions, one finally obtains the present day contribution of these GWs to the critical density, 
i.e. dOgw/dlnA; = k/(3Hq) dp gW;0 /dln k 



G-cnd 



dlnfc 625 V a ) I H 



kj(k)l(k) , (35) 



recalling the definitions introduced before: k = k/{a cn ^H en d)> q = <z/(a C nd-ffcnd) and introducing the dimensionless 
quantities 

J{k) = l6H 2 nd x^k 6 J(k) = H* nd P 2 J{k) , (36) 

rt».x r+i , N (l-u 2 ) 2 a 3 

T(k) = / dg / d M P c (g)P c ( vV + fc 2 - 2fcg/i) — -i ' \ ... , (37) 



i 



3/2 ' 

if + fc 2 - 2&g> ) 



where, in the argument of the power spectrum, q and k 
must be considered as function of q and k according to 
the simple relations introduced above. 

Expression Eq. (l35l may be further simplified when 
assuming a standard thermal history, i.e. radiation dom- 
ination between the reheating epoch and the epoch of 
matter radiation-equality shortly before recombination. 
In this case 

£^"» r '^(iilv)"" ! 

xkj(k)i n (k), (38) 

where 51 7i o — 5 x 10~ 5 is the present density param- 
eter in radiation and where we have defined In(fc) by 
In(fc) = 10 22 X(k) since T{k) is quadratic in Vq- When 
reheating occurs long after the end of inflation, i.e. for 
Xf s» x T \i ^> Xini ~ %end which implies a r h 3> a on d, the 
expression in the curly brackets of Eq. (131 1) approaches 
asymptotically x^ h . In this case ^ approaches unity. 

Finally, it may also be interesting to have an expression 
giving the associated present day GW frequency. One 



finds 

1 /3 

/o - 5xlo5Hz Kl^v) ' (39) 

In Fig. [2j we show the result of a numerical integra- 
tion of the above formula for two different reheating tem- 
peratures T r h = 10 9 GeV and T r h = 10 7 GeV, respec- 
tively, along with projected sensitivities of next gener- 
ation gravitational wave experiments, namely advanced 
LIGO THill, BBO and DECIGO ||. The prediction 
for the above gravitational wave signal is shown by the 
thick solid line. At this stage, it is important to stress 
that the above calculation is carried out in the linear ap- 
proximation and cannot describe the non-linear aspects. 
In particular, there exists a range of modes [k n i, fc max ] for 
which growth is sufficiently effective that these modes be- 
come non-linear by the time reheating starts. This range 
of wave numbers is discussed in the next Section. Here, 
we wish to stress that the linear growth Sk <x a, as used 
in the above calculation, is not adequate anymore in the 
non-linear regime. This has impact on the above calcu- 
lation in the range [k n \, fc max ]. In order to show how the 
calculation could be affected, we have plotted in thick 



7 




10° 10 1 10 2 10 3 10 4 10 5 10 6 10" 2 10° 10 2 10 4 10 6 

frequency [Hz] frequency [Hz] 



FIG. 2: Energy density in gravitational waves (in unit of the present day critical density) as a function of frequency. Left 
panel: T r h = 10 9 GeV; right panel: T r h = 10 7 GeV. The thick solid line shows the result of the linear calculation of the GW 
signal. The thick short dashed line shows the same calculation, but when stopped at the onset of the non-linear stage (see 
text). Various non- linear signal (discussed in Section are shown: the dotted line ("collapse") shows the signal associated 
with the collapse of halos that have gone non-linear by reheating; the upper limit "merger" presents an estimate of the possible 
amplification of the collapse signal due to tidal interactions; the dashed-dotted line ( "turbulence" ) shows the signal expected 
from the turbulence produced by the evaporation of non-linear structures at reheating; and the star symbol ("evaporation") 
shows the signal produced by the evaporation of the collapsed halos at reheating. The long dashed line annotated "inflation" 
shows the signal expected from the direct generation of GWs in m 2 cj> 2 inflation. On the left hand side, the various lines indicate 
the projected sensitivities of next generation gravitational wave detectors. The horizontal dotted line gives the limit imposed 
by big-bang nucleosynthesis. 



dashed line in Fig. [5] the result obtained when one shuts 
off the emission of gravitational waves when mode k or 
|fc — q\ become non-linear. This curve, which sets a strict 
lower limit to the gravitational wave signal departs from 
the above linear calculation in the range [fc n i , fc max ] ■ 

In Fig. [51 we also plot several estimates of gravitational 
wave production in the non-linear regime, to be discussed 
in Section IIVI We also plot the signal expected from the 
direct production of gravitational waves in m 2 2 inflation 
(represented as the long dashed line annotated "infla- 
tion"}. We have calculated this signal following the study 
of [2J]. The change of spectral evolution (from nearly 
scale invariant at small frequencies to oc k~ 2 ) occurs at a 
frequency fc m i n /(27r), with k min — a r h-ff r h- Indeed, modes 
with k > fc m in have re-entered the Hubble radius during 
inflaton domination in the pre-reheating epoch and suf- 
fered redshifting while those modes with k < k m i n have 
re-entered the horizon in the radiation dominated era (up 
to the modes with very small k that re-enter in the matter 
dominated era after matter-radiation equality, of course). 

It is clear that the magnitude of the linear signal is 
weak for the reheating temperatures considered here. It 
is important to recall, however, that such reheating tem- 
peratures lie in the upper range of constraints set by the 
impact of gravitinos on big-bang nucleosynthesis and that 
in the absence of a definite underlying particle physics 
model in the inflaton sector, there is no particular scale 
for T rn . The reheating temperature could be much lower, 
in which case the strength of the signal would increase 



(roughly as T rh , see above) and the frequency would 

1 /3 

decrease (as T rh ), coming closer to the range of maxi- 
mum sensitivity of future instruments. 

IV. GENERATION OF GRAVITATIONAL 
WAVES IN THE NON-LINEAR REGIME 

A. Collapse of halos 

After perturbations have grown in the linear regime 
to Spk/p ~ 1 they separate from the Hubble flow and 
collapse to form halos. In order to calculate the grav- 
itational wave signal coming from this phase, we as- 
sume that the halos collapse and virialize instantaneously 
when the fluctuation on the scale considered exceeds the 
threshold 8 C ~ 1.7. With this approximation, a halo of 
mass M emits gravitational waves at a well-defined time 
and at a single pulsation 2tt/, which is estimated as the 
reciprocal of the dynamical timescale of collapse i C oii- 

The energy density emitted during the collapse of 
perturbations corresponding to mass between M and 
M + AM is 

j coll _ i coll^coll, { Q>cou\ , An \ 

dPg W ^ dn h £ gw tcoul-^— I • (40) 

In this expression, is the luminosity emitted by one 
halo of mass M at collapse, dn™ 11 represents the number 
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density of halos of mass between M and M + dM at the 
time of collapse, and t co \\ ~ is the collapse 

timescale (with p n the density of the halo at the time of 
collapse, see below). The factor (a co n/a now ) 4 accounts for 
the redshifting of the signal through cosmic expansion. 
The value of a co ii as a function of wave number will be 
determined further below. A halo of mass M corresponds 
to a structure characterized by a comoving wave-number 
k such that (assuming spherical symmetry) 



M 



end 



Therefore, the above equation can be rewritten as 



dfl 



coll 



drC 11 dlnM u 



din/ 3H£ dlnM din/ 



sw ''coll 



a 



(41) 



(42) 



with / the comoving frequency of emission of gravita- 
tional waves, / ~ (27rf co ii) -1 , to be evaluated further 
below. 

The luminosity C^ 1 can be evaluated as follows [25| . 
The typical amplitude of a gravitational wave at point x 
(the origin of the coordinates being chosen at the location 
of the emitting object) is given by the quadrupole formula 
which reads 



G 
~2 



(43) 



where n% = Xj/|x|. The quantity in brackets is the trace 
free part of the quadrupole tensor J^ . Since spherical 
motions do not generate gravitational waves, we need to 
assume that the halos are not spherically symmetric. A 
rough estimate of the quantity Zy is 



h 



2 / PhViVj 
hi 



2GM 
R 



(44) 



where R is the typical size of the halo (physical radius) at 
virialization and the virialized velocity u v ir — \J 2GM j R. 
The luminosity of the source is thus given by 



"•coll 



G 



G M° 



(45) 



where we have written h ~ h/T with T ~ R/v V i V . The 
radius R is related to the wavenumber of the fluctuation 
by R = 2na coU /k. 

Formally, the number density of collapsed halos at a 
given time can be derived from a Press-Schechter mass 
function. With the above approximation at that all halos 
of mass M collapse instantaneously when the rms mass 
fluctuation on that scale exceeds 5 C , one derives the num- 
ber density of halos of mass between M and M + dM at 
collapse as 



dn 



coll ^ p(Qcou) 
M 



d In M 



&cnd 



M 



dlnM 



(46) 



with p(a co ii) the mean density of the Universe at the time 
of collapse. This relation is nothing but the differen- 
tial version of the equation rih/V = p/M (V is the vol- 
ume) which expresses the fact that, according to Press- 
Schechter, at a given time, most of the matter is con- 
tained in the halos that collapse at that time. 

The rms mass fluctuation om is related to au , the root- 
mean-square of the Fourier amplitude of density pertur- 
bations via 



M 



d 3 fc' 



(27T) 



(47) 



with k the comoving wave number of the perturbation 
enclosing a mass M defined in Eq. (j4T|) . In the last equa- 
tion an approximation of simple k-spacc filtering has been 
applied, i.e. only modes k' < k contribute to the average. 
Using Eq. (|4"7|) one finds 



cm 



(48) 



Collapse of mass scale M occurs when um — S c , hence 
at a collapse scale factor a co n defined by 



Qcoll 
Q-cnd 



= 55 c 7' c (fc)" 1/2 r 



(49) 



At collapse, the typical over-density <5 co ii ~ 150 (in the 
spherical collapse model). The present day frequency is 
related to the comoving wave number of the perturbation 
as 



/n 



1 a coll „ /■, , x sl/2 

-Hcoll (1 + Ocoll) , 



2tt ao 
3.8 x 10 3 HzA: 



1/3 

Vl0 9 GeVy V ID 



1/4 



r^rvj^i (50, 



Combining all the above, one ends up with 



dn 



coll 

gw 



din/ 



6 x IQ^^o—V 5 ^ 



2 x 10 _23 fc- 2 



T,, 



10 9 GeV 



4/3 



io- 11 



5/4 

(51) 



Maximum emission is provided by the largest scales that 
go non-linear at reheating, i.e those with wave number 
A: n i such that a co ii(fe n i) = a rn - According to Eq. (|4T)|) . this 
corresponds to 



knl 



' 10 9 GeV/ llO- 1 



-1/4 



(52) 



where we approximated V^(k n \) = 10 11 in accord with 
Fig. [T] The frequency at maximum emission is thus 



f£° n ~ 3.7 x 10 2 Hz 



10 9 GeV 



(53) 



This falls in an interesting range considering the LIGO 
or BBO gravitational wave detectors. Nevertheless, the 
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signal is fairly weak: at the frequency of peak emission, 
one finds dfi^jj'/dln/ ~ I0~ 20 independent of reheating 
temperature, a signal which would be only very hard to 
detect. The signal as a function of frequency is repre- 
sented in Fig. [2] in dotted line, annotated "collapse". 

B. After initial collapse of halos 

We have so far assumed that efficient GW emission 
in the non-linear regime happens only during the col- 
lapse phase of halos, up to the point where virialization 
is completed. It is not clear if post-virialization rota- 
tion/vibration of halos, potentially enhanced by occa- 
sional tidal forces acting on the halos due to passing 
by of other halos, may lead to further efficient emis- 
sion of GWs. The question of how much more GWs are 
emitted after halo formation may only be addressed by 
complete numerical simulations and is beyond the scope 
of this paper. However, we may attempt to estimate 
it, introducing an efficiency parameter e < 1, for emis- 
sion of GWs after virialization. We may formally relate 
the generated GWs at times between the end of collapse 
and reheating Pg W co11 to p^ 11 in Eq. (|4T))) by making the 
following replacements in Eq. (|30|): p™ 11 ' -> Pg W co11 ' , 
dn h (M) -> dn> coll (A/) = dn h {M){a co \\/af to take into 
account the fact that the signal is emitted after the col- 
lapse, £g^ n — > £gw° n = e £ g w' as discussed above, and, 
finally, t co u -> £> co n = £ C oii(a/a col i) 3 / 2 , (a col i/a ) 4 -> 
(a/ao) 4 , again to take into account that the emission time 
is now different. Inserting this into Eq. (|4TJ)) yields 

d^^dp^^y \, (54) 

implying a potentially large enhancement of the signal 
when a r h 3> a co ii, provided e is not too small. Using 
this in conjunction with Eq. (|51|l one finds that for non- 
negligible e the signal may be dominated by the emis- 
sion of GWs at T r h from the first structures which had 
formed. One would then expect a signal df2 gw /dln/ ~ 
2 x 10- 18 e(T rh /I0 9 GeV)- 2 ('P c /10- 11 ) 5 / 2 with a rather 
strong dependence on the reheating temperature, which 
would make it potentially detectable (provided e is not 
too small, of course): for T r h ~ 10 6 GeV, the signal be- 
comes of order 10" 12 in the lCT 2 Hz range. Assuming 
indeed that the typical pulsation corresponds to the Hub- 
ble scale at the time of emission, one finds a frequency 
/ >co11 - 30Hz(T rh /10 9 GeV) for the peak emission at 
T r h- This signal is indicated by the "merger" upper limit 
in Fig. [5J given the uncertainty on the value of e. 

C. Evaporation of halos 

Around the epoch of reheating, the self-gravitating ha- 
los evaporate due to emission of radiation by the decaying 
inflaton. Due to the smallness of the inflaton-radiation 



coupling, the radiation does not scatter on the inflaton 
particles contained in the halos, but escapes freely with- 
out communicating its pressure to the inflaton halo. In 
this picture, each halo behaves as a source of radiation 
with a typical emission timescale = T^ 1 corresponding 
to the inflaton lifetime. The mean free path of radiation 
against scattering with itself is A ~ l/(a 2 T), with a a 
typical coupling constant and T the temperature of the 
radiation, therefore A <C and the radiation can be 
approximated as an instantaneously thermalized ultra- 
relativistic fluid. The radiation emitted by the evaporat- 
ing halo expands and accelerates under its own internal 
pressure until it interacts with the winds emitted by sur- 
rounding halos. We will address this phase in the fol- 
lowing subsection and for the time being, we estimate 
the amount of gravitational waves produced during the 
evaporation of the halos, taken individually. 

Since the radiation escapes on a timescale t csc = R/c, 
with R the radius of the virialized halo, and since 
R -C H^ 1 , the hierarchy i csc -C remains satisfied 
all throughout evaporation. This implies that the energy 
of radiation contained in the halo at any time, Mi osc /r0 
is always smaller than M, so that its gravitational influ- 
ence can be neglected. 

As the halo evaporates, its virial radius increases in 
inverse proportion to M. In order to see this, one first 
notices that the dynamical timescale of evolution of the 
halo, idyn = { K Ph/3)~ 1 remains much smaller than the 
Hubble time (rep/3) -1 / 2 during the evaporation process, 
since ph p. Therefore the structure adjusts itself on 
a timescale which is much smaller than the evaporation 
timescale = (H~^ the Hubble time at reheat- 

ing). In other words, the halo undergoes a series of quasi- 
static equilibria as the mass decreases, the virial relation 
Mv 2 ~ 2GM 2 /R being verified at each step. Through 
adiabatic expansion, the rms velocity decreases as l/R, 
hence one can check that dM/M — —dR/R. 

During evaporation, the luminosity can be expressed 
as [H[ 

C^~eG(^pj . (55) 

Here e < 1 is an efficiency factor which represents a mea- 
sure of the asphericity of the radiation wind, analogous 
to the factor introduced in Section TlVBI Using the pre- 
vious arguments, the amount of energy density emitted 
in gravitational waves on an evaporation timescale ap- 
proximately reads: 

G fdM\ 2 , 
dp gw ~ £ j(-^-) T <t> dn h \ A {M) 

~ e^M 2 H lh dn hllh (M) (56) 

where drihirh is the number density of halos of mass com- 
prised between M and M + dM at the time of reheat- 
ing and 7J r h — 1/(2^) denotes the Hubble scale at re- 
heating. It is important to note that all halos evaporate 



10 



on the same timescale r^, therefore they all contribute 
to the gravitational wave signal at the same comoving 
pulsation ~ a T bT$ ■ Due to the quadratic dependence 
on mass scale M the signal is dominated by the largest 
structures evaporating. This is given by the typical mass 
scale which collapses shortly before reheating. It may be 
obtained by employing Eq. with a co \\ « a r h yielding 
the final simple result 



gw 



£-O 7 , P C (fcnl) 3/4 



(57) 



independent of reheating temperature. The typical fre- 
quency of emission corresponds to the (redshifted) Hub- 
ble scale of reheating (divided by 2n) and is given as 
before by / >co11 - 30 Hz (T rh /10 9 GeV), in an interesting 
range from the point of view of detection. 

The above signal may be quite substantial, in partic- 
ular, for P^ ~ 10~ n as implied by Fig. Q] one obtains 
n gw ~4x 10~ 13 e. Of course, the signal may be smaller 
than our estimate in case velocities are below the speed 
of light and/or £<1, 



D. Dissipation after reheating 

As explained above, once the radiation is emitted, it 
will shock against the radiation generated by surround- 
ing halos. All in all, this injects a large amount of energy 
on all scales < L s = 2ira I \ 1 /k n i, where k n i denotes as be- 
fore the smallest comoving wave number of perturbations 
that go non-linear by reheating. On general grounds, one 
expects that this injection of energy stirs a turbulent cas- 
cade on scales < L Sl which decays through damping at 
small scales. This cascade is characterized by a power 
spectrum of kinetic energy in modes (eddies) of various 
wave numbers, which can emit gravitational radiation in 
the very same way that growing perturbations do during 
the linear phase. This has been discussed in detail by 
Kosowsky et. al [10| and we simply adapt their calcu- 
lation to the present framework. We may assume here 
that a fraction of order unity of the total energy density 
at reheating is injected in the cascade, i.e. kinetic energy 
of the radiation fluid is of order the rest energy. The 
stirring scale L s can be evaluated by using Eq. (j^)) , with 
«coii = irh, which yields 



2P 1/4 P" 1 

Z ' C rh 



(58) 



The stirring scale is much smaller than the timescale ~ 
H~L over which energy is injected into the plasma, so 
that one may consider the turbulence to be stationary 
over this timescale t^. Following [l(|, one then finds 
that the gravitational wave signal generated is: 



dfL 



din/ 



2 x 10 



-1 / T \3/r\ -7/2 
-7 I T <f> \ I ^ s \ I J 



rh 



2 x 10 



io- 



H Th ) \f s 
£ 



(59) 



independent of reheating temperature. The typical fre- 
quency is defined as 



f 1 fl rh / r 2\-V3 



10 d Hz 



I 10 9 GeV/ I 10 11 



-1/6 



(60) 



Although not as strong as that resulting from the evapo- 
ration of halos, this signal remains substantial and the 
typical frequency can easily shift down to the inter- 
esting range for detection if the reheating temperature 
T rh < 10 7 GeV. This signal is indicated in Fig. [2] by the 
dashed-dotted line annotated "turbulence" . 



V. CONCLUSION 

In this paper we have calculated the gravitational wave 
signal associated to the growth of sub-horizon perturba- 
tions between the end of cosmic inflation and the begin- 
ning of a radiation dominated early Universe in m 2 <f> 
chaotic inflation. The growth of metric fluctuations for- 
mally arises from a preheating like instability and can be 
interpreted, in the range of wave numbers of interest, as 
the gravitational instability of a pressure-less fluid in a 
matter dominated era. The over-density Spk / p of these 
perturbations grows linearly with scale factor and reach 
non-linearity to re-collapse and form inflaton halos after 
only moderate expansion of the Universe. 

Though these small-scale inflaton halos later evapo- 
rate when the inflaton decays at reheating, they may 
nevertheless lead to the emission of gravitational waves 
in an interesting frequency range luq ~ 10~ 5 — 10 8 Hz 
for current ground-based and planned satellite gravita- 
tional wave detectors, such as Advanced LIGO, DECIGO 
and BBO. We have therefore analyzed the gravitational 
wave emission in detail during five distinct phases: (a) 
growth of the linear perturbations, (b) first collapse of 
halos when reaching non-linearity, (c) subsequent non- 
linear evolution of halos, (d) evaporation of halos due to 
inflaton decay, and (e) cosmic turbulence during the first 
epochs of radiation domination, respectively. We note 
here that phase (a) has recently also been studied by 
Ref. 2l| . An exact calculation of (a) indicates that gravi- 
tational wave emission is weak, peaking at values of order 
tt gw ~ 10" 22 T r ; 4/3 ,ifthc magnitude of curvature pertur- 
bations on small scales as predicted in chaotic inflation, 
i.e. Pq ~ 10 -11 , is assumed. The estimate of gravita- 
tional wave emission during phase (b) is similarly small. 
Promising gravitational wave signals of il gw > 10~ 12 for 
~ 10~ n , could however, result during phases (c) (d) 
and (e) , depending on the value of the reheating temper- 
ature. Generally speaking, the lower the reheating tem- 
perature, the longer the phase of growth of fluctuations 
and the larger the band of wave numbers turning non- 
linear, and consequently, the larger the signal of gravi- 
tational wave emission. The exact value of these signals 
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depends on the details of the structure formation process phericity, etc. A reliable estimate may therefore only be 
including tidal disruption of substructures, rotation, as- given when full numerical simulations are performed. 
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